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Abstract 

x; 

We present the general expressions for the resummation, up to next-to-leading loga- 
rithmic accuracy, of Sudakov-type logarithms in processes with an arbirtrary number of 
hard-scattering partons. These results document the formulae used by the authors in several 
previous phenomenological studies. The resummation formulae presented here, which are 
valid for phase-space factorizable observables, determine the resummation correction in a 
process-independent fashion. All process dependence is encoded in the colour and flavour 
structure of the leading order and virtual one-loop amplitudes, and in Sudakov weights as- 
sociated to the cross section kinematics. We explicitly illustrate the application to the case 
of Drell-Yan and prompt-photon production. 
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The perturbative QCD calculations of a large class of infrared and collinear safe observables 
are sensitive to Sudakov effects. Some classical examples of these observables are the e + e~ energy- 
energy correlation in the back-to-back region pQ, the cross section for Drell-Yan production of 
lepton pairs in hadron collisions j2], and several e + e~ hadronic event shapes in nearly two-jet 
configurations [3]. 

The Sudakov effects appear when the observable is defined and/or measured close to the 
exclusive boundary of its phase space. We generically denote by y (y > 0) the kinematical variable 
that measures the distance from the exclusive boundary, so that the Sudakov region is specified by 
)/< 1. When Sudakov sensitive observable are computed as power series expansions in the QCD 
coupling as, the perturbative series involves terms of the type a"gL k (k < 2n), where L = —Iny. 
These double logarithmic terms are due to final-state radiation of soft and collinear partons, and 
are a distinctive feature of any short-distance dynamics that is governed by an underlying gauge 
field theory. Since L ^> 1, the presence of logarithmically-enhanced terms spoils the convergence 
of the fixed-order expansion in as{Q 2 ), even if the observable is controlled by a typical hard- 
scattering scale Q whose value is large (such that as{Q 2 ) < 1). The predictivity of perturbative 
QCD can be recovered by reorganizing the perturbative series according to the degree of divergence 
of the various logarithmic terms, and then by performing a systematic resummation, to all orders 
in as, of the contributions that are leading- logarithmic (LL), next-to- leading logarithmic (NLL), 
and so forth. 

Resummed calculations up to NLL accuracy are available for several production cross sections 
in hadron collisions (see the list of references in Sect. 5 of Ref. [3]), and for many hadronic event 
shapes in e + e~ annihilation (see e.g. Refs. [HIE]) and in deep-inelastic lepton-hadron scattering 
(see e.g. Refs. [U]). The inclusion of resummed Sudakov effects increases the theoretical accuracy 
of perturbative calculations, by extending their applicability to wider phase-space regions and 
reducing the uncertainty coming from yet uncalculated higher-order terms. This brings about 
relevant improvements in phenomenological applications, as shown by the studies carried out 
in recent years [T. For example, in e + e~ annihilation the use of resummed calculations has 
become the standard procedure in the comparison with data on hadronic event shapes J8J: these 
calculations allow one to extend the perturbative treatment towards the two-jet region where 
statistics is higher; they also allow investigations of hadronic physics at the interface between 
perturbative and non-perturbative phenomena. In hadron collisions, resummed calculations often 
lead to a considerable reduction in the scale dependence of the perturbative predictions, as in the 
case of top quark production at the Tevatron and bottom quark production at HERA B j§] EH E] • 

In recent years, different groups (KLOS [T2], BCMN [3 CH|, BSZ [13]) have been working to 
develop resummation formalisms that are process- independent and observable- independent. The 
aim is to obtain generalized resummation formulae that depend on universal coefficients, and that 
are applicable to different hard-scattering processes and different classes of observables within 
the same process in terms of a minimal amount of information on the specific observable to be 
computed. We have explicitly checked that our generalized resummation formulae (which are 
presented here) reproduce known NLL results for several quantities, such as the thrust [El El and 
C-parameter [15J distributions in e + e~ annihilation, the cross sections for the production of lepton 
pairs, vector bosons j2] and Higgs bosons [TH] in hadron collisions, the structure functions [T71 [TH] 
in deep-inelastic lepton-hadron scattering at large values of the Bjorken variable. We used this 
formalism to derive the NLL resummed results of Ref. for the production of heavy quarks 
and prompt photons in hadron collisions. However, a general description of the formalism has 
never appeared in the literature. The purpose of this work is to fill this gap. Here we only give a 
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brief illustration of our generalized resummation formulae. More details on the formalism and its 
derivation are given in a forthcoming paper. 

The paper is organized as follows. We first consider QCD hard-scattering processes without 
hadrons in the initial state. We discuss the kinematic properties of the observables to which our 
resummation formalism applies. Then, we present our generalized resummation formula up to 
NLL accuracy. The explicit formula is expressed in terms of factorized final-state factors (Jj) 
and interference terms (A^ mt '- ) ). Then, we discuss the more general case of hard scattering in 
hadron collisions and in processes with tagged hadrons in the final state. Here the corresponding 
resummation formulae include additional initial- and final-state factors (Aj). We briefly illustrate 
the application of the general formalism by sketching the derivation of the resummation formulae 
presented in Ref. ^U]. Finally, we summarize our main results. 

We begin our presentation by considering a generic infrared- and collinear-safe cross section 
er (or a related observable) in a hard process that does not involve hadrons in the initial state 
(for instance, hadron production in lepton collisions or in heavy-boson decays). We suppose that 
the calculation of o at the leading order (LO) in QCD perturbation theory involves m final-state 
QCD partons with four-momenta {p{\ = p±, . . . ,p m . For simplicity of presentation, we also limit 
ourselves to considering the case of massless {pf = 0) QCD partons (quarks, antiquarks and 
gluons). Using a shorthand notation, we write the LO contribution a^ L °^ to the cross section as 

a (LO) = J d$((T;{p . }) \M^ \{ Pl })\\ (1) 

where M^ L °^ is the corresponding LO matrix element, and \M <yLO ' > \ 2 denotes the squared matrix 
element summed over the colours and spins of the final-state QCD partons. The kinematics of the 
cross section are fully described by the phase-space factor d<£>(<7; {pi})- It includes the phase-space 
contributions for the production of the final-state particles as well as any additional kinematics 
information (definition of jets, event shapes, energy flows, ...) that is necessary to precisely 
define the cross section a that we want to evaluate. The phase-space dependence on a is briefly 
indicated by the notation d$(<j). In particular, d$(<r) depends on the generic kinematic variable 
y that controls the distance from the Sudakov region. We assume that the LO term 
well-behaved* (not singular) as y — > 0, while higher-order terms contain logarithmically-enhanced 
contributions of relative order a^L 2n . The dependence of a on the momenta of non QCD partons 
(7, Z°, W ± , H, . . .) is always understood. 

Note that the Sudakov logarithms in a do not necessarily occur by approaching the true phys- 
ical phase-space boundary. These logarithms can also appear inside the phase space of certain 
observables* . Indeed, logarithmically-enhanced terms may arise jTH] also if the phase-space bound- 
ary for a certain number of partons lies inside that for a larger number, or if the observable itself 
is defined in a non-smooth way at some perturbative orders. In these cases, a^ 10 ^ in Eq. (Q) has 
to be regarded as the lowest-order contribution at which those partonic boundaries appear. 

The practical feasibility of performing the resummation of the Sudakov logarithms at all pertur- 
bative orders depends on the capability of properly approximating the higher-order contributions 
to a. The approximation regards both the QCD dynamics (i.e. the matrix elements) and the 

Un practice, we consider the case in which all the LO invariants PiPj are of the order of the hard scale Q 2 when 

$A notable example in e + e~ annihilation is the C-parameter distribution, which has Sudakov logarithms in the 
vicinity of C = 3/4 [T^J- Other examples are discussed, for instance, in Refs. |rJll2U). 
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cross section kinematics. As for dynamics, since the Sudakov limit singles out multiple radia- 
tion of soft and collinear partons, we can exploit the universal (process-independent) factorization 
properties of the QCD multiparton matrix elements in the infrared (soft and collinear) region (see 
e.g. Refs. [21] and [221 )• As for kinematics, we restrict our study to a (large) class of observables, 
whose phase space is factorizable. By phase-space factorization we precisely mean the following. 
At higher perturbative orders, we consider the contribution to a from the final-state radiation 
of additional partons with momenta c?i, ■ ■ ■ , m the Sudakov limit y — > 0, these momenta are 
kinematically forced to become soft or collinear to the momenta {pi}, and the cross section o is 
called factorizable if the corresponding phase space d$(er;pi, . . . ,p m , q l7 . . . , q k ) behaves as 

k 

d$(a;p 1 ,...,p m ,q 1 ,...,q k ) —> d$(a;{pi}) [dq] u(a, {p { }; qj) , (2) 

y j=i 

where d$(a; {pi}) is the LO phase space, [dg] = Yij d 4 <?j<5+(<?|)/(27r) 3 is the phase-space contribu- 
tion from the unconstrained emission of the additional partons^ with on-shell momenta qi, . . . ,q k , 
and on the right-hand side we have neglected relative corrections that vanish in the soft and 
collinear limit. The function u(a, {pi}; qj) is called Sudakov weight. It depends on the kinematical 
definition of the cross section a (such dependence implicitly embodies the dependence on y), on 
the LO parton momenta {p^ and on a single (soft and collinear) final-state momentum qj. The 
right-hand side of Eq. (j2J) implies that the kinematics dependence on the soft and collinear mo- 
menta is fully factorized: it is factorized with respect to the LO phase space and, moreover, there 
are no correlations between those momenta, since each Sudakov- weight factor depends on a single 
momentum qj. 

Note that the momenta {p{\ on the right-hand side of Eq. (J2J) are not precisely the momenta of 
the LO partons on the left-hand side. The former exactly coincide with the latter in the soft limit 
qj — >• 0. When some of the momenta qj are not soft but collinear to the momentum of one of the 
LO partons, say the parton i, the momentum pi on the right-hand side is obtained by reabsorbing 
the longitudinal-momentum recoil produced by the collinear radiation. 

As a consequence of the infrared and collineary safety of a, the Sudakov weight fulfils the 
following important property: 

u(a, {pi}; q) = 1 when q = , or q = (1 — z)pi for i = 1, . . . , m . (3) 

Note that infrared and collinear safe observables are not necessarily factorizable. A classical 
example of non-factorizable observables are jet rates when the jets are defined by the JADE 
jet-finder algorithm |53] . Moreover, phase-space factorization is typically not achievable in the 
space of the kinematic variables where the cross section is originally defined. To overcome non- 
factorizable phase-space constraints, it is often necessary to introduce a conjugate space. For 
instance, the constraints of energy or transverse-momentum conservation are usually factorized 
by respectively performing Mellin (or Laplace) or Fourier transformations, and by working in the 
TV-moment [21 El 121] or impact-parameter PQ space. Thus Eq. ([2]) can be valid either in the original 
space or in a properly defined conjugate space. In the following, y generically stands for either 
the original Sudakov variable or the variable conjugate to it (more precisely, the inverse of this 
conjugate variable) in the conjugate space. 

§We are treating the partons as distinguishable particles. It the partons j = 1, . . . , k were identical, [dq] should 
be multiplied by a Bose-symmetry factor of 1/kl. 
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To proceed further, we require one additional kinematics property on the observable to be 
resummed. In the Sudakov limit, the energy flow accompanying the LO hard scattering has to 
be suppressed uniformly with respect to its emission direction. To be precise, we consider the 
behaviour of the Sudakov weight u(q) (from now on, u(a, {pi}] q) is briefly denoted by u(q)) when 
y — > at fixed value of q. We write the on-shell four- momentum g M = q) in terms of its energy 
qo = u) and a three-dimensional vector q of unit lenght (q 2 = 1), whose components parametrize 
the emission angle. When y — > 0, we thus require that u(q) — > and, without loss of generality, 
we can always assume, to dominant logarithmic accuracy, that u(q) is approximable by a step 
function: 

u(q) ~ 6 (u; max - u) , w max (?/; {p^}, q) ^ — ^ , (4) 

where, the upper bound u> max on the radiated energy depends on the momenta {pi}, on the emission 
angle q and on the Sudakov variable y. In the Sudakov limit, we then require 

dln^ max (?/; {pi},q) 

amy y — > 

where the power A is positive, independent of y and {p{\ and, in particular, independent of the 
radiation angle q. Equations (JH) and (jSJ) state in a formal way that, in the Sudakov limit, the 
parametric suppression rate of the energy flow emitted from the LO partons has to be uniform 
with respect to the radiation angle. 

Note that, by requiring the property in Eqs. (J3J) and (j5J), we exclude from our resummation 
treatment observables such as away-from-jet energy flows, and, in general, the so-called non-global 
observables 

The general kinematic properties of phase-space factorization (see Eq. and uniform sup- 
pression of the energy flow (see Eqs. (@J) and are sufficient to obtain our generalized resumma- 
tion formula. The all-order cross section is generically denoted by o u . Here o u can be either the 
original cross section a, or the corresponding cross section in the conjugate space where Eq. (J2J) 
applies (in this case, a is eventually computed by perfoming the inverse transformation of o~ u to 
the original space). The cross section a u is written as 

a u = a^ + at ] , (6) 

where the resummed component a^ es ^ contains all the Sudakov logarithms, while CT^ is well 
behaved (finite or vanishing) order by order in as when y — > 0. Thus, can reliably be 

evaluated by truncating its perturbative expansion at the first few perturbative orders. In practice, 
af n) can be obtained from the fixed-order computation of a u by subtraction of the terms already 
included in cr^ res ^ at the same fixed order. 

The resummed component is given by 

^ rcs) = / d<Ma; fe}) \M^°\{ Pl })\ 2 , (7) 

where d$ u (a; {pi}) is either the LO phase space d$(cr;{pj}), or its version in the conjugate 
space. The expression (J7J) is completely analogous to the LO expression in Eq. (0). Sudakov 
resummation is simply achieved starting from the LO result and performing the replacement 

^To be precise, we allow Eqs. (@J and © to be violated in angular regions of vanishing solid angle (for instance, 
when q is exactly parallel to a LO momentum pi). 
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|j^(LO)|2 _^ |j\/f( L °)| 2 £( M ). The generalized effective form factor embodies the dependence 
on the Sudakov logarithms to all perturbative orders. Since the Sudakov limit y — > can for- 
mally be regarded as the limit u — * 0, the presence of logarithmically-enhanced terms in S(it) is 
identified by contributions that order by order in Qg are formally divergent when u(q) — > 0. 

The Sudakov logarithms exponentiate, and £ has the following structure 

£ ~ C(a s )exp{g(a s ,L)} = [1 + a s Ci + . . .] exp {L g 1 (a s L) + g 2 {a s L) + . . .} . (8) 

The coefficient factor C(as) is independent of u (or y) and is due to hard virtual radiation. Its 
perturbative coefficients C\,C%,... depend on the specific cross section a, but, since they are 
not logarithmically enhanced, they can be computed process by process at some finite perturba- 
tive orders. The exponent Q(ag,L) contains the logarithmically-enhanced terms. The function 
Lgi(ctsL) resums the LL contributions a$L n+1 in the exponent, the function g 2 (asL) resums all 
the NLL contributions agL n , and so forth". 

The explicit NLL resummation formulae we are going to present have the structure of Eq. (jSJ), 
with the only difference that the functions C(ois) and Q(as,L) are matrices in the flavour and 
colour indices of the hard-scattering partons, so that exponentiation has to be understood in 
formal sense. The expression of the form factor £ up to NLL accuracy is 

The first factor on the right-hand side of Eq. Q contains all the LL terms and part of the NLL 
terms. It is given by the product of the jet functions Jj(w), and there is a jet function Jj for each 
final-state parton i in the corresponding LO process. The function Ji(u), which generalizes the 
jet function of Refs. [21 EH EH], embodies all the logarithmic terms produced by multiple radiation 
of partons that are collinear (either soft or not) to the direction of the momentum pi of the 
LO parton i. The factorization in single-parton factors, Jj, is a consequence of the independent 
character (which follows from colour coherence) of QCD collinear radiation. 

The remaining factor on the right-hand side of Eq. Q contains NLL terms and subleading 
logarithmic contributions. The radiative factor A^ mt ^(w) generalizes the analogous radiative factor 
introduced in the case of prompt-photon hadroproduction [TQ]. It embodies all the quantum- 
interference effects produced by non-collinear (large-angle) soft-gluon radiation. In particular, 
this factor is sensitive to the colour correlations due to the colour flow dynamics of the LO hard 
scattering. 

The function Ji(u) has the following resummed expression: 

(2tt) 3 ""^ ' ^ p { q 

where the functions Pi(z), which depend on the flavour (i = q, q, g) of the parton i, are related to 
the Altarelli-Parisi splitting functions 



lnJ 4 (u) =4vr J 5 + {q 2 ) {u{q) - 1) A {a s {2{\ - z iq ) Pi q)) Pfa) , (10) 



i + z 2 „ r 2 , . i 



f 



z 2 + (l-z) 2 



'ID 



P q (z) = P- q (z) = Cfj-j , P g (z) = C A - 2 + z(l - 

"in our definition of LL, NLL, etc. terms, we are referring to the logarithmic hierarchy of the various contri- 
butions to the exponent Q(as,L) (i.e. to InE) in Eq. JSJ. Our systematic resummation procedure thus differs 
from the ones that refer to the expansion of E in successive logarithmic towers, such as dgL 2n , a 1 gL 2n ~ 1 , etc.. In 
particular, our LL and NLL contributions include more logarithmic terms than those in the first two logarithmic 
towers of E (see, for instance, the discussion in Sect. 5 of Ref. |2fi|). 
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and A(as) is the QCD coupling as denned in the bremsstrahlung scheme [T7], and is the related 
to the MS coupling as by the NLO relation 



A(a s ) = a s 




K = CA\V---\--N f . (12) 



The 'energy fraction' Zi q is defined with respect to a four-momentum nf as 

ma 

z iq =l-^-. 13 

riiPi 

Indeed, to obtain our resummed formulae, we have introduced m auxiliary momenta nf. These 
auxiliary momenta are arbitrary, with the only constraint of being time-like (nf > 0) and hard** . 
The LL terms in Ji{u) do not depend on the definition of these auxiliary momenta. The NLL 
dependence of Ji(u) on rii is cancelled by the dependence on rij in the remaining factor on the 
right-hand side of Eq. Q, so that our expression for £(«) is independent of rij up to corrections 
that are beyond the NLL accuracy of the present formalism. The main motivation for introducing 
the auxiliary momenta is to factorize collinear radiation in the jet functions Ji(u) without 
the introduction of explicit angular boundaries. In practical calculations, the definition of the 
momenta can be adjusted to simplify the evaluation of the integral in Eq. (fTU|) . 

The NLL contribution A^ mt ) to the form factor is given by 

A (int.)( M ) = v( M ) V(u) , (14) 

where V(w) and V(m) are matrices acting onto the colour indices of the LO partons. The explicit 
expression of V(-u) is 

V(.)-P. e X p(g/ d ^Ml_£M) £^^1 , ( 15 , 
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where is the renormalization scale, to be chosen of the order of the hard-scattering scale Q. The 
sum runs over the labels i, j = 1, . . . , m of the LO partons and T£, (c = 1, . . . , N% — 1) 
are the corresponding colour charges: T? is the colour matrix^ in the fundamental (adjoint) 
representation if the parton i is a quark (gluon). The operator P z denotes z-ordering in the formal 
expansion of the exponential matrix, and V(w) is obtained from Eq. ()15j) by simply replacing P z 
with P z , the ordering operator that acts in the opposite order. The notation (u(q)) in Eq. (fTo'j) 
stands for a properly defined average of the Sudakov weight u(q). Parametrizing the light-like 
four-vector as q^ = u;(l,q), the average is performed over the angular directions q at fixed 
value q = u = (1 — z)^,r of its energy q . In practice, exploiting Eqs. (JU) and ©, and neglecting 
terms beyond NLL accuracy, we simply have 

(«(?)> = © {y x - (i - *)) • (is) 

In Eq. ( |9*|) . the NLL colour matrix A^ 1111 -* acts onto the colour indices of the m-parton matrix 
element Mn{{pi})- This hard matrix element is independent of u and is perturbatively computable 
as a power series in as(^ 2 R ): 

M H ({ Pl }) = M(°\{ Pi }) + as^ R ) Mg\{pi}) + 0(a 2 s ) , (17) 



**By hard, we mean that the invariants nf and riipj are of the order of the hard scale Q 2 when y 
ft The colour charges are defined according to the notation in Sect. 3.2 of Ref. |27j . 
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where is the LO matrix element and is the hard part of its one-loop virtual corrections. 
Since M# does not containt logarithmically-enhanced terms, it can be evaluated by truncation of 
its perturbative expansion at some fixed order. In Eq. @, the dependence of M H on the colour 
indices is represented by the colour vectors \M H ) and (M H \, which are defined according to the 
notation in Sect. 3.2 of Ref. |27j . 

The algebraic complications due to the non-trivial colour structure of the NLL contributions 
are straightforwardly overcomed when the number m of LO partons is m = 2 or m = 3. In these 
cases the colour algebra can be carried out in closed form, since the colour matrices Y^cTfTj in 
Eq. (|15j) are simply proportional to the unity matrix in colour space. The proportionality relations 
are (see Appendix A in Ref. |2Tj ) 

£7\ C T 2 C = -Ci = -C 2 (m = 2) , (18) 

c 

£7\ C T 2 C = i(C 8 -C 1 -C7 2 ) (m = 3) , (19) 

c ^ 

where Cj is the Casimir invariant of the parton % (Cj = Ca if the parton % is a gluon, C, = Cp if 
the parton i is a quark or an antiquark). When m — 3, Dc^~2^~3 and 53 c Tg are obtained from 
Eq. (|T9"j) by permutation of the parton indices {1,2,3}. In the case of processes with m > 4 LO 
partons, the colour algebra cannot be handled without additional information on the hard matrix 
element Mfj{{Pi}), since the colour can flow in many different ways through the hard scattering. 
In general, the evaluation of X^T^TJ (and of Eq. (JISI)) requires the diagonalization of a linear 
combination of mim — 3) /2 independent colour matrices (see Appendix A in Ref. [2*7]). 

The NLL resummed calculations of Ref. ^2] deal with the non-trivial colour structure in 
the hard-scattering of m = 4 LO partons by considering 'soft anomalous dimensions' of gauge- 
dependent Wilson line operators. The Wilson line operators are introduced and properly defined 
on a pro cess- dependent basis. In this respect, the integrand in the exponent of Eq. ([To]) can be 
regarded as a universal (process-independent and observable-independent) soft anomalous dimen- 
sion matrix (in colour space), T, of soft non-collinear gluons radiated by hard scattering of an 
arbitrary number (m > 4) of partons. The explicit expression of T, 

a s T({ Pt ,n t })=as £ E?T?7 ' (2 ° } 

i^j c \Ptyj) 'H 11 j 

is gauge independent, though dependent on the auxiliary vectors rij. The Altarelli-Parisi splitting 
function Pi(z) in Eq. (jlOJ) controls the resummation of the Sudakov logarithms produced by 
collinear (soft and hard) radiation, and, by analogy, T controls the resummation of the Sudakov 
logarithms produced by soft non-collinear radiation. Obviously, the definition of the boundary 
between the collinear and non-collinear regions is quite arbitrary. In Eq. (J2U|) , this arbitrariness is 
somehow parametrized by the dependence on the auxiliary vectors nj. For instance, using colour 
conservation, Y^i=\Ti — |2~7j . Eq. (}2T)j) can be rewritten as 

a S T({ Ph n i })=asEET^ln^-2astc i ln^ , (21) 

ijtj c KPiPj) i=\ lb i h^i 

where /ij are abitrary scales (e.g. fii = fiji). The second term on the right-hand side is proportional 
to the unity matrix in colour space, and therefore it can be absorbed in a corresponding redefinition 
of the jet functions Jj in Eq. (JHJ). 
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We now consider the general case of cross sections in hadron collision processes. In these 
processes, the hadronic cross section a is obtained by convoluting the partonic cross sections 
cr ai,a 2 with the parton densities f ai /hi( x i, H F ) and fa 2 /h 2 ( x 2,^ 2 F) of the colliding hadrons with 
momenta Pi and P 2 . 

<?= I dx l dx 2 fai/hifa, fip) f a2 / h2 (x 2 ,H 2 F ) O axm , (22) 

11,12 

where [ip is the factorization scale (to be chosen of the order of the hard-scattering scale Q), and 
the sum Yl ai ,a 2 runs over ^ ne flavours (ai,a 2 = q,q,g) of the incoming partons with momenta 
Pi = XiP\ and p 2 = x 2 P 2 . 

At the LO in QCD perturbation theory, the partonic cross section <J ai ,a 2 has the same structure 
as in Eq. (0), with m LO partons i (i — 1, 2 are the incoming partons and i = 3, . . . , m are the 
outgoing partons). The Sudakov logarithms at higher orders can be produced by soft and collinear 
radiation emitted from either the outgoing partons or the incoming partons in the LO hard 
scattering. In general, to kinematically factorize the Sudakov effects in the parton densities from 
those in the partonic cross section, it is necessary to consider the iV-moments (Mellin moments) 
fa/h^i/J'p) °f the parton densities, 

/a/fc,iv(/4) = J dx X N ~ l f a /h(x, fJ.%) , (23) 

and the corresponding iV-moments cr^ 1 ^ 2 of the partonic cross section. 

In the following we limit ourselves** to considering the cases in which hard radiation collinear 
to the incoming partons is suppressed in the Sudakov limit, so that the Sudakov logarithms are 
produced by soft radiation (at any angles) and hard radiation collinear to the outgoing partons. 
This simplifies the presentation of the resummed formulae, since the Sudakov effects do not change 
the flavours ai, a 2 of the incoming partons. 

Our NLL resummed formulae apply to the iV-moments c^ 1 ^ 2 of the partonic cross sections 
that fulfil the factorization property of Eq. (j2J) in a properly defined conjugate space. We still 
require the properties in Eqs. (j3J) and (jSJ). Infrared and collinear safety implies Eq. (J3J) in the case 
of radiation collinear to the outgoing partons i — 3, . . . , m. When i — 1, 2 is an incoming parton, 
Eq. (j3J) is modified by a kinematical rescaling factor, which simply takes into account that we are 
considering the iVj-moments of the partonic cross section. We have 

u(q) = z^" 1 ~ exp{-(iVi - 1)(1 - z)} when q = (1 - z)p { for i = 1, 2 , (24) 

where the approximate equality is valid in the soft region (1 — z 1) we are interested in. The 
all-order partonic cross section has the same structure as in Eq. and its resummed component 
is obtained as in Eq. (JJJ). The only difference is that the (final-state) form factor has to be 
replaced by a more general radiative factor A aii iY 1 ;a2,Af 2 (' u )- The latter is given by an expression 
similar to Eq. ©, apart from two simple modifications that regard the LL and NLL terms, 
respectively. 

The modification of S(tt) at the LL level is obtained by performing the replacement 



II Mu) — A auNl (u) A a2 , N2 (u) ]J Mu) ■ (25) 




ttWe thus exclude observables such as, for example, the QT-distribution of Drell-Yan lepton pairs, where powers 
of IiiQt are produced also by radiation of hard quarks and gluons that are collinear to the colliding partons. 
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In other words, A aijNl . a2:N2 (u) is obtained from by supplementing the product of the final- 
state jet functions Ji(u) with an initial-state factor A aij Ar(w) for each incoming parton i — 1,2. 
The initial-state Sudakov factor A^(m), which generalizes the analogous iV-moment factor of 
Refs. [21 El US] , resums the Sudakov logarithms produced by soft partons emitted collinearly to 
the LO incoming parton i = 1,2. As in the case of the parton density f a /h,N (/^f); the radiative 
factor A Ni {u) depends on the factorization scale \ip and on the factorization scheme used to define 
the partonic cross section. The NLL resummed expression of A^-(u) in the MS factorization 
scheme is 

lnA aiM (u) = 47r/^5 + (<7 2 )M<?)-^)) A(a s (2(l - z iq ) Pi q)) P t (z iq ) 

The term on the right-hand side of the first line of Eq. (|2*Bj) is completely analogous to the 
right-hand side of Eq. (JTUJ) apart from the replacement (u(q) — 1) — > (u(q) — Ui(z iq )), where 
Ui(z) = u(q = (1 — z)pi). The Sudakov limit typically forces the parton distributions towards the 
large- N (large-x) region, and the term in Eq. (J26|) matches the sensitivity of the parton distribution 
fa i /h,N i {.^ 2 F) to large logarithms, IniVj, of the Mellin index 

The modification of S(-u) at the NLL level regards the interference term A^ int ^ or, more 
precisely, its components V(-u) and V(w) in Eq. (fTlj) . Equation (fTTjj) gives V(m) as an exponential 
of the soft anomalous dimension matrix T in Eq. (J2*U|) . When going from processes with no 
initial-state partons to processes in hadron-hadron collisions, we have to take into account that 
the parton momenta pi and pi have to be crossed from the final to the initial state. As for 
the anomalous dimension T, the crossing simply amounts to the following analytic continuation, 
ln(piPj) — > ln(piPj)+iir, of the terms with i = 1 or i = 2 and j > 3. Performing such a replacement, 
and using colour-charge conservation, jyjLs^i = ~ (^T + ^2% we obtain the overall replacement 
to be applied to T: 

r(fe,^})^r({ Pi ,^})±4m^(T 1 c + T 2 c )(T 1 c + T 2 c ) . (27) 

c 

Here, the signs + and — regard the replacements to be applied in the evaluation of V(ii) and 
V(it), respectively. Note that the substitution on the right-hand side of Eq. (|2Tj) is effective only 
in the case of hadron-hadron collisions with m > 4 hard-scattering partons at LO. In fact, when 
m = 2 we have + T% = 0, and when m = 3 we have S c (Tf + T 2 C )(T 1 C + T^) = Cf so that the 
substitution in Eq. (j2*7|) changes V(it) by a pure phase factor that is cancelled by its complex 
conjugate phase factor in V(w). 

The replacement in Eq. ()27j) , to be applied to the soft anomalous dimensions in hadron-hadron 
collisions, has a direct and interesting physical interpretation. As is well know from QED (see e.g. 
Ref. [28 ), Coulomb-type virtual exchanges produce infrared- divergent Coulomb phases that affect 
any scattering amplitudes. Being them phases, they cancel in the evaluation of inclusive cross 
sections. The QCD analogue [22] of the QED Coulomb phases are non-abelian Coulomb 'phases', 
which are colour matrices. They lead to non-abelian infrared divergences that cancel, as in QED, 
in infrared- and collinear-safe observables with incoming massless partons [2II1I22J- However, in 
the non-abelian case, the cancellation mechanism of the infrared divergences can produce residual 
(and non-trivial) finite contributions. The Sudakov logarithms produced by the i7r-term on the 
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right-hand side of Eq. (|27|) are a manifestation of these residual* finite effects. Very soft (and, 
hence, infrared divergent) non-abelian Coulomb-type interactions do cancel in a. On the contrary, 
non-abelian Coulomb-type interactions of gluons that are harder (and, hence, infrared finite) than 
bremsstrahlung gluons do not cancel, since they are trapped by the colour fluctuations produced 
by the radiated bremsstrahlung gluons. 

We note that our resummation formulae apply also to processes in which the partonic cross 
section has to be convoluted with the partonic fragmentation functions of hadrons that are tagged 
in the final state. The resummation formulae for these processes are simply obtained by multi- 
plicatively introducing a factor of A au Ni( u ) f° r eacri final-state parton i whose momentum has to 
be convoluted with the Ni moment of the corresponding fragmentation function. The fragmenta- 
tion factor A a . j ]sf i (u) is the same as for parton densities (i.e. Eq. (|26Jl ). provided the fragmentation 
functions are defined in the MS factorization scheme. 

To illustrate the use of our generalized resummation formulae, we briefly sketch their appli- 
cation to the hadroproduction of Drell-Yan lepton pairs [2] and of prompt photons |iQ: • In both 
cases, the hadronic cross section is obtained by the factorization formula in Eq. (|22p. and v^S* 
(S = (Pi + P2) 2 ) denotes the centre-of-mass energy. The hard-scattering scale is the invariant 
mass Q of the lepton pair in the first case, and the transverse energy E? of the photon in the 
second case. The variable y that parametrizes the distance from the Sudakov region is respectively 
given by y — 1 — Q 2 / S and y = 1 — 4E%,/S. We are interested in the corresponding inclusive total 
(i.e. integrated over the rapidity of the observed final state) cross sections in the Sudakov limit 
j/ < 1. We thus consider the N- moments of the cross section in Eq. (|22|) . The N- moments are 
defined with respect to (1 — y) at fixed values of Q and Et, respectively. This sets N\ = N2 = N 
in the corresponding partonic cross sections cr^ 1 ^ 2 . Since the moment variable N is conjugate to 
y through the Mellin transformation, the Sudakov limit corresponds to 1/N — > (i.e. N — > 00). 
It is not difficult to prove that the factorization property of Eq. (J2J) is valid in iV-space for both 
processes. 

In the Drell-Yan process, the LO hard-scattering subprocess is q(pi) + qij>2) — > ££'(Q): the 
number of LO partons is m = 2 and the LO kinematics set Q 2 = 1p\P2- The Sudakov weight is 



which fulfils Eqs. (jlj and (J5J) with A = 1. Since the LO process has no final-state partons, the 
Sudakov resummation factor Ady,n is the product of A^ mt ^ and two initial-state factors, A qii N 
and Aq 2i 7v- To apply our generalized resummation formulae, it is convenient to choose the auxiliary 
vectors as n\ = = p\ +P2- We thus have A^ int -) = 1 (see Eqs. (fTlj) and (fT5)l) and A qi: N = Aq 2 ^ 
(see Eq. (|26Jl ). Moreover, since u(q) = Ui(z iq ), the term on the right-hand side of the first line of 
Eq. (|26|) vanishes. Therefore, from the second line of Eq. (J26j) we finally obtain the known NLL 
result for the Drell-Yan process 0. 

In the prompt-photon process (see the second paper in Ref. PHI); there are two independent 
LO hard-scattering subprocesses, ax(pi) + ^2(^2) — ¥ 03(^3) +7(^7); with {ai = q,a 2 = q,a 3 = g} 
and {at = q, (12 = g, 03 = q}- The number of LO partons is m = 3 and, in the Sudakov limit the 

*The i7r-term in Eq. (|27|l is not the absolute overall effect of Coulomb- type interactions in hadron-hadron 
collisions. It is the relative effect produced by the Coulomb-phase mismatch between processes with no initial-states 
hadrons and processes in hadron-hadron collisions. Thus, it has conveniently been introduced by the corresponding 
analytic continuation procedure. 




(28) 



10 



LO kinematics set 2E^ — 2pxp 3 ~ 2p 2 ps — P\Pi- The Sudakov weight is 

% (g)=exp(-iV^) , (29) 

which fulfils Eqs. (@J) and © with A = 1. The Sudakov resummation factor A aia3 _^ a37)J v is the 
product of two initial-state factors (A 01) jy and A a2i jv), one final-state factor (J a3 ) and the NLL 
factor A^ nt '). To explicitly evaluate these factors, it is convenient to choose the following auxiliary 
vectors: ri\ — p\ + p$, n 2 = P2 + P3, n 3 = p\ + p 2 + p^- Using this choice, it is straightforward to 
check that the term on the right-hand side of the first line of Eq. (|2l)j) gives subleading (beyond the 
NLL accuracy) contributions. As for the colour algebra, we can simply use Eq. (|19j) . Performing 
the phase-space integrals in Eqs. (jTUjl . (fT5)l and using Eq. (J2EJ), we obtain the explicit resummed 
results anticipated in Sect. 4.2 of the the second paper in Ref. jTU] and implemented in the 
phenomenological study of Ref. [3U] . 

We have discussed a generalized formalism to perform the resummation of Sudakov logarithms 
in QCD hard-scattering processes. We have presented explicit resummation formulae up to NLL 
accuracy. The formulae are observable- independent and process- independent. The dependence 
on the observable is universally encoded in the one-particle Sudakov weight u(a, {p^}; q). The 
dependence on the process is completely specified by flavour, colour charge and kinematics of the 
LO partons, the LO matrix element and its hard virtual corrections at one-loop order. This is the 
minimal amount of pro cess- dependent information that is necessary in any calculations at fixed 
perturbative order. Within a specific process, the formalism is applicable to a large class of QCD 
observables that are specified by some kinematic properties, such as phase-space factorizability. 
Phase-space factorization has already been exploited in the literature to perfom resummation of 
several observables in processes controlled by LO hard scattering of two and three QCD partons. 
The extension of the NLL resummation techniques to multiparton processes requires a formalism 
able to deal with the radiation pattern of non-collinear soft gluons emitted by the colour flow 
dynamics of the underlying hard-scattering. Available NLL resummed calculations of some spe- 
cific cross sections in four-parton processes treat the colour flow dynamics on a process-dependent 
basis. As for processes with higher number of LO partons, no NLL resummed calculation has 
been presented so far. Our NLL formalism applies to arbitrary processes with any number of 
hard-scattering partons and with arbitrary colour flow dynamics. This opens prospects of phe- 
nomenological applications to multijet events at present and future high-energy colliders. 
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